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ABSTRACT 

The radio source Sgr A* at the center of our Galaxy is believed to be a 2.6 x 10 6 Af© black hole which 
accretes gas from the winds of nearby stars. We show that limits on the X-ray and infrared emission 
from the Galactic Center provide an upper limit of ~ 8 x 1O" 5 M yr" 1 on the mass accretion rate 
in Sgr A*. The advection-dominated accretion flow (ADAF) model favors a rate <, 1O" 5 M yr" 1 . In 
comparison, the Bondi accretion rate onto Sgr A*, estimated using the observed spatial distribution of 
mass losing stars and assuming non- interacting stellar winds, is ~ 3 x 10 Af yr . There is thus rough 
agreement between the Bondi, the ADAF, and the X-ray inferred accretion rates for Sgr A*. We discuss 
uncertainties in these estimates, emphasizing the importance of upcoming observations by the Chandra 
X-ray observatory (CXO) for tightening the X-ray derived limits. 

Subject Headings: accretion, accretion disks - black holes - Galaxy: center - radio sources: Sgr A* 
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1. INTRODUCTION 

Stellar kinematics show that there are « 2.6 x 1O 6 M 
within w 0.015 pc of the Galactic Center (Eckart & Genzel 
1997, Ghez et al. 1998), centered on the radio source Sgr 
A* (Menten et al. 1997). The most plausible explanation 
is that Sgr A* is a « 2.6 x 10 M accreting black hole. 

Sgr A* is believed to accrete the winds from nearby 
(~ 0.1 pc) massive stars (Krabbe et al. 1991). Hydro- 
dynamical simulations of the Galactic Center region, as- 
suming 10 randomly distributed point sources (to mimic 
the nearby stars) and a 10 6 M© black hole, yield accretion 
rates w 1 - 2 x 1O" 4 M yr" 1 (Coker & Melia 1997; here- 
after CM) . If we naively increase this by a (Bondi capture) 
factor of 2.6 2 to account for the actual mass of the central 
object, the theoretically predicted accretion rate becomes 
- lO" a M yr" 1 (but see §3). 

The bolometric luminosity of Sgr A* is observed to be 
<, 10 37 ergs s" 1 (Genzel et al. 1994, Narayan et al. 1998). 
For an accretion rate of ~ 1O _3 M yr -1 , this corresponds 
to a radiative efficiency of ~ 10" 7 ! A possible explana- 
tion for the low luminosity of Sgr A* is that the accretion 
occurs via an advection-dominated accretion flow (ADAF; 
Narayan, Yi & Mahadevan 1995, Manmoto et al. 1997, 
Narayan et al. 1998; see, e.g., Melia 1992, 1994, Mas- 
tichiadis & Ozernoy 1994, Falcke 1996, Beckert & Duschl 
1997 for alternative models of Sgr A*). All ADAF mod- 
els in the literature, however, require an accretion rate 
of < 1O" 5 M yr" 1 (see Table 2 of Quataert & Narayan 
1999, which gives accretion rates in Eddington units; their 
estimates need to be multiplied by a factor of 0.055 for 
Mq yr" 1 ). If CM's black hole mass scaled Bondi capture 
estimate of ~ 1O _3 M yr" 1 is taken at face value, there is 
a discrepancy of a factor of <; 10 2 in mass accretion rates. 

In §2, we sharpen this discrepancy by using spectral 
observations of the Galactic Center to constrain the accre- 
tion rate of gas at large distances from the central black 
hole. The argument is nearly independent of the accre- 
tion model employed. We then calculate Bondi capture 



estimates of accretion onto Sgr A*, explicitly using the 
observed spatial distribution and wind properties of stars 
in the Galactic Center (§3). In §4 & §5 we conclude and 
assess the agreement between various theoretical estimates 
of the accretion rate in Sgr A* . 

2. OBSERVATIONAL CONSTRAINTS ON M 

We use X-ray and IR limits on the emission from the 
Galactic Center to derive constraints on the accretion rate 
of optically thin or thick gas onto Sgr A*. We take the 
distance to the Galactic Center to be 8.0 kpc. In what 
follows, 77i 5 denotes the mass of the black hole in units of 
2.6 x 1O 6 M , R denotes the physical radius in the accre- 
tion flow, r (= R/7.38 x 10 11 ms cm) denotes the radius 
in Schwarzschild units, and rh is the accretion rate in Ed- 
dington units, where M e dd = O.O55m,gM yr" 1 . 

2.1. X-ray Constraints on Optically Thin Accretion 

If the gas accreting onto Sgr A* is optically thin, it will 
emit X-rays by bremsstrahlung emission. The observed 
limits on X-ray emission from Sgr A* can therefore be 
used to constrain the accretion rate. The accretion rate 
of the gas is related to the density, p, by the continuity 
equation 



M = AirR 2 pv fS ri H ri v = M a (r/r ) p . 



(1) 



where r Q ~ 10 5 , the "outer" radius of the flow, is roughly 
the Bondi capture radius of the black hole. The radial 
velocity of the gas is related to the free fall velocity, Vff , 
by \v\ = rjyVff and the vertical scale height of the gas is 
H = t/hR- For a pure Bondi flow, T) v — t]h = 1, while for 
an ADAF, r] v ~ a (the viscosity parameter) and r/H ~ 1/2. 

Equation (|T|) allows for the possibility that the mass 
accretion rate may decrease with radius due to an out- 
flow/wind (in which case the assumed geometry is a 
"thick" equatorial inflow with a bipolar outflowing wind; 
cf Blandford & Begelman 1999). The accretion rate at the 
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outer edge of the flow (i.e., the mass supply rate) is M Q 
and the parameter p measures the strength of the outflow. 

We parameterize the electron temperature in the accre- 
tion flow by T e = ri T T {r/r )- a , where T = 10 12 /r o is ~ 
the proton virial temperature at r Q and a (~ 1 at large 
radii) determines the radial profile of the electron temper- 
ature; t\t ~ 1 for both Bondi accretion and an ADAF. 

The bremsstrahlung luminosity at frequency v is given 

by 



Am-jH I e v RrdR, 



where e v = 2.6 x 10 iu /9 2 T e _1 ^ 2 exp[— hv/kT e ] ergss^ 1 cm -3 
is the bremsstrahlung emissivity (e.g., Rybicki & Lightman 
1979; we have ignored a weak frequency-dependent Gaunt 
factor). Carrying out the integral in equation (||) using 
the above expressions for the density and temperature of 
the gas, and approximating the exponential cutoff in e„ by 
a step function, we find (in ergs s _1 Hz -1 ) 



(2) 
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where m = M a /M e dd and 
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(4) 



Observations at frequency v are dominated by emission 
from the radius r v . In particular, if the X-ray energies of 
interest are less than the minimum electron thermal en- 
ergy in the flow (which is obtained at r a and is ~ 1 keV 
for r Q ~ 10 5 ), the bremsstrahlung luminosity is dominated 
by emission from the largest radii in the flow. In this case 
r v = r Q and the parameters p and a, which determine 
the radial profiles of the density and electron tempera- 
ture, are irrelevant; they enter only in the normalization 
of the emission, not in the exponent of any parameters of 
the problem. As a result, observations of soft X-rays pro- 
vide effective constraints on M D , the rate at which matter 
is supplied to the accretion flow on the outside. They can- 
not, however, constrain the flow structure closer to the 
black hole (e.g., the nature or strength of an outflow). 

Parameterizing the observational constraints on emis- 
sion from the Galactic Center via vL v 55 10 34 L34_ergs s _1 
at hv = v\ keV [y = 2.4 x 10 17 z^i Hz) equation (||) can be 
inverted to give a constraint on the accretion rate of the 



gas (taking 2p 
simplicity) : 



a/2 ss 1 in the denominator of eq. 13] for 
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where r 5 — r o /10 5 . 

The parameters r/ v , 7Jh , and i]t must all be < 1 . 
thermore, the Bondi capture radius in Sgr A* for 

1 Genzel et al. (1997) report a possible detection of Sgr A* 
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material moving at ~ 10 3 km s _1 is ~ 10 5 Schwarzschild 
radii so that the accretion flow must originate at r 5 ~ 1. 

ROSAT observations of the Galactic Center indicate 
that the X-ray luminosity of Sgr A* at « 1 keV is <; 10 34 
ergs s _1 , in which case L34 pa V\ w 1 and r v pa r a 
(but, see Narayan et al. 1998 for a discussion of the 
effect of uncertainties in the absorbing column; we take 



N, 



H 



6 x 10 22 cm 



By equation (|^) this implies 
M ^ 8 x 1O _5 M0 yr -1 ; this is a strict upper limit on 
the accretion rate of optically thin gas in Sgr A* . 

ASCA observations of the Galactic Center (Koyama et 
al. 1996) give a 2 - 10 keV limit of 10 35 ergs s -1 . Taking 
v\ w 5, equation (||) gives a limit of M v < 10 _4 Af Q yr -1 
at r v « 0.2r o ~ 10 4 , where M„ = M (r„/r ) p is the ac- 

SIGMA ob- 

35 



-1/0 1/0 



cretion rate at the radius r v w 

servations (Vargas et al. 1998) give vL y ^ 2 x 10 35 ergs 
s _1 at vi k, 75, which requires M v < lO _4 Af0 yr -1 at 
r v m 0.01r o ~ 10 3 . If there arc no winds emanating from 
Sgr A* at r ~ r (i.e., p«0) then the ASCA and SIGMA 
observations require M ^ lO _4 Af0 yr -1 . This may be an 
even stronger limit than that derived from the ROSAT ob- 
servations, since it is insensitive to the absorbing column. 

2.2. Infrared Constraints on Optically Thick Accretion 

The preceding subsection assumes that the gas is opti- 
cally thin at large radii. We believe that this is the most 
likely scenario. It is worth, however, examining obser- 
vational constraints imposed on the accretion rate of an 
optically thick disk. 

For a standard Shakura-Sunyaev disk (e.g., Frank, King, 
& Raine 1992), the luminosity is given by 



vL,. 



16Tr 2 h cos(i)v 4 



RdR 



exp(hu/kT(R)) -V 



(6) 



where i is the inclination of the disk, R ou t and Ri n are the 
inner and outer radii of the disk, and T(R), the effective 
temperature of the disk, depends on the accretion rate (cf 
Frank et al. 1992). 

Menten et al. (1997) have obtained an upper limit of 
ps 10 35 ergs s _1 on the 2.2 fim emission from the Galactic 
Center. This emission strongly constrains the properties 
of any geometrically thin disk. It is straightforward to 
numerically integrate equation (^) and obtain the 2.2 /1 
emission from the disk. Requiring this to be Ss 10 35 ergs 
s~ 4 yields a constraint on the parameters Af and ri n . For 
i - 60°, taking Af w 1O" 3 M yr" 1 implies r m k, 10 4 . If 
M w 10 _5 Mq yr -1 the constraint is slightly less stringent, 
r.i n <; 3 x 10 3 (note that, to satisfy all existing IR and ra- 
dio limits on emission from Sgr A* typically requires even 
larger values of r^). 1 

This implies that, for any reasonable inclination and 
accretion rate, a thin disk cannot extend much inside 
r.i n ~ 10 4 without violating the IR limits on emission 
from the Galactic center. Even if we assume that from 
the Bondi capture radius to ri n the accretion is via a ge- 
ometrically thin disk, and then for r < ri n it is optically 
thin, our limits on Af from the previous subsection are 



at 2.2 (im, at a flux level above Menten et al.'s upper limit (indicating possible 
variability in the source). This need not represent a detection of an optically thick disk, but could be either (1) singly Compton scattered 
synchrotron emission, (2) synchrotron emission from non-thermal electrons, or (3) optically thick accretion disk emission. 7/ interpreted as the 
latter, the upper limits in this paragraph become (approximate) equalities. 
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not changed substantially. This is because the X-ray con- 
straints on the accretion rate are quite insensitive to the 
precise outer boundary of the optically thin region (cf. eq. 
H where the limit is cx r^ 1 ^ 4 ). 

3. BONDI CAPTURE ESTIMATES 

There are ra 10 He I stars (Krabbe et al. 1991) pro- 
jected within 10 arcsec of the position of Sgr A* (Menten 
et al. 1997). The winds from these stars result in mass loss 
rates ranging from about 8 x 10 -5 to 8 x 10~ 4 M© yr _1 
and might provide most of the matter accreted by Sgr A* . 

Using the distribution of the He I stars projected on 
the plane of the sky, their mass loss rates, M w , and their 
wind speeds, V w (Najarro et al. 1997), we have estimated 
the mass accretion rate onto a 2.6 x 1O 6 M black hole. 
The line-of-sight (z-coordinate) offset of each star relative 
to Sgr A* (the black hole) is unknown, and we assigned 
a z-distance randomly to each star, assuming a Gaussian 
distribution. The stellar wind material passing within a 
radius Ra — 2GMbh/V% of Sgr A* was assumed to be 
accreted. This resulted in a contribution to the accretion 
rate for a given star, M acci given by 

Mace = \M W - y/l-Ryr^j , (7) 

where r is the three-dimensional distance of the star from 
Sgr A* (Note that all stars in the Najarro et al. sample 
have r > Ra)- 

Using equation (0), a star such as IRS 16C with M w ra 
8 x 10~ 5 M & yr _1 , a wind speed of 650 km s _1 , and a true 
distance from Sgr A* of ra 2 arcsec (its projected distance 
is 1.4 arcsec) would contribute about 1.4 x 10 -5 Mq yr" 1 
to the total mass accretion rate of Sgr A*. On the other 
hand, the most powerful He I star near Sgr A*, IRS 13E1, 
which loses mass at a rate of 8 x 10 -4 M© yr _1 with a wind 
speed of 1000 km s -1 , has a projected distance from Sgr 
A* of 3.7 arcsec and thus makes a maximum contribution 
to the mass accretion rate of only 5 x 10~ 6 Mq yr . 

For each trial distribution of stars, obtained for differ- 
ent realizations of the random z-coordinate distances from 
Sgr A* , we calculated the total mass accretion rate as the 
sum of the accretion rates for each star. This assumes no 
wind-wind interactions. Varying the standard deviation of 
the z-coordinate distribution, er z , we found average total 
mass accretion rates of M a = 5, 3, and 2 x 10~ 5 M yr -1 
for a z of 2, 5, and 8 arcseconds, respectively. 2 If the He 
I stars have a <r z ra 5, as suggested by the projected spa- 
tial distribution of stars near Sgr A* (Genzel et al. 1997; 
Ghcz et al. 1998), then the total mass accretion rate for 
non-interacting stellar winds onto a 2.6 x 10 6 Mq black 
hole should be M ra 3 x 10" 5 M yr" 1 . 

This accretion rate is noticeably smaller than the value 
of ~ 10~ 3 M© yr -1 obtained by scaling CM's numerical 
results from a 10 6 M Q black hole to a 2.6 x 10 6 M Q black 
hole (cf. §1). The main reason for the discrepancy seems to 
be that the spatial distribution of point sources they take 
(e.g., ±3Ra) is a reasonable approximation of the observed 
stellar distribution only for a 2.6 x 10 6 -M© black hole, not 
a 10 6 M Q hole. That is, for v w ra 700 km s -1 and a black 
hole mass of 1O 6 M , 3Ra is only ra 1", which is smaller 

2 For an "average" wind speed of 700^7 km s — 1 , this corresponds to 



than the typical projected distance of the mass losing stars 
from Sgr A* (on the other hand, for a 2.6 x 1O 6 M black 
hole, 3Ra corresponds to ra 3"). One should therefore 
not scale CM's result up by a factor of 2.6 2 . This leaves a 
smaller residual discrepancy, perhaps due to wind-wind in- 
teractions, which they treat in their hydrodynamical sim- 
ulations, but which we neglect. 

4. DISCUSSION 

In §2 of this paper, we showed that ROSAT X-ray ob- 
servations of Sgr A* provide a firm limit on the accretion 

rate at large radii, namely M a ^8x X0~ 5 r) v L^ Mq yr -1 
at r ra 10 5 , where r/ v is the ratio of the radial velocity of 
the gas to the free fall velocity and L34 is the soft X-ray 
(ra 1 keV) luminosity of the source. We obtain similar 
limits at r ~ 10 4 and r - 10 3 using ASCA and SIGMA 
data, respectively. It is important to emphasize that these 
limits correspond to the X-ray radiation being dominated 
by bremsstrahlung emission. If part of the X-rays are pro- 
duced by Comptonization, as in some models, then the 
accretion rates must be even lower than the limits derived 
here. 

The X-ray luminosity of 10 34 ergs -1 used in the above 
scalings is itself an upper limit. There is considerable dif- 
fuse emission from the Galactic Center (Koyama et al. 
1996), and the luminosity of the accretion flow is uncer- 
tain. Another uncertainty is the absorbing column, which 
could lie anywhere from the value used here, 6 x 10 22 cm~ 2 
(Watson et al. 1981), to 1.5 x 10 23 cm- 2 (Predehl & 
Trumper 1994). If the column is closer to the higher value, 
then the unabsorbed luminosity of Sgr A* could be larger 
than we have assumed. High angular resolution observa- 
tions with CXO (formerly AXAF) should solve both prob- 
lems. 

Sgr A* probably accretes most of its mass from the 
winds of nearby massive stars (Krabbe et al. 1991). We 
have estimated the Bondi accretion rate onto Sgr A* using 
the observed spatial distribution of these stars, along with 
their estimated mass loss rates and outflow velocities (§3). 
The resulting accretion rate is M Q ra 3 x 1O _5 M yr -1 
for plausible values of the unknown spatial distribution of 
stars along the line of sight. This estimate of M Q should 
perhaps be interpreted as a lower limit, as there may be 
other sources of matter for the accretion flow, and wind- 
wind interactions and post-shock cooling may increase the 
fraction of the wind material accreted. 

The upper limit on M Q from the X-ray observations and 
the lower limit on M Q from the capture of stellar winds 
are roughly compatible, requiring 77^, ~ 1 for the accre- 
tion flow and L34 ~ 1 for the (unknown) X-ray luminos- 
ity. If the accretion onto Sgr A* occurs viscously via an 
ADAF, then rj v ra a. For values of a ra 0.01 seen in 
numerical simulations of thin disks (Hawley, Gammie, & 
Balbus, 1996), there is a substantial discrepancy between 
the Bondi and X-ray inferred accretion rates, while for the 
values of a ra 0.1 — 0.3 favored by Narayan et al. (1998; see 
also Esin et al. 1997), the two estimates are comparable. 

Blandford & Begelman (1999) have recently argued that 
only a small fraction of the mass supplied to an ADAF 
reaches the central black hole, with most of it being lost 
to an outflow/ wind. It is important to emphasize that this 

a z of 2, 5, & 8 2 R A . 
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alone will not modify the conclusions of this paper. Our 
primary constraint on the accretion rate utilizes the ~ 1 
keV X-ray flux, which is dominated by bremsstrahlung 
emission from large radii in the flow (r « r ~ 10 5 ); at 
such radii, any wind from the accretion flow has little ef- 
fect. 

Gruzinov (1999) has argued that radial heat conduction 
would significantly reduce the mass accretion rate in Bondi 
flows. This explanation is also unlikely to effect the results 
of this paper. The density at ~ ro is a boundary condition 
set by the stellar winds. Gruzinov's flow would therefore 
produce the same X-ray luminosity from radii r ~ ro as the 
models we have considered (at large radii ~ rn, Gruzinov's 
reduced accretion rate is due to a smaller radial velocity, 
not a lower density). 

Finally, it is worth emphasizing that our paper does not 
explicitly address the reason for the low radiative efficiency 
of Sgr A*. This is determined by the physical conditions 
close to the black hole (r <, 10 — 100), while our analysis 
focuses on larger radii. 

5. CONCLUSION 

With current data, we believe that there is rough agree- 
ment between the Bondi accretion rate (as estimated in 
this paper), the limits inferred from X-ray and IR obser- 
vations, and the value (or, more precisely, the upper limit) 
favored by ADAF models. If this straightforward interpre- 
tation is correct, CXO should detect an X-ray flux compa- 
rable to that of ROSAT, despite its significantly improved 
angular resolution. On the other hand, if CXO substan- 
tially lowers the soft X-ray flux from the Galactic Center 
there would be an inconsistency between the Bondi rate 
and the values favored by the X-ray observations. Given 
the uncertainties in the accretion rate estimates, it is worth 
spelling out what we believe to be the two plausible, if 
somewhat mundane, explanations if such an inconsistency 
is discovered. 



The first assumes that the Bondi accretion rate esti- 
mates are correct and M ^ 3x 10~ 5 A/ Q yr -1 . This could 
be reconciled with the limits inferred from a significantly 
smaller soft X-ray flux only if Njj is larger, by a factor of 
few, than the value we have assumed (6 x 10 22 cm -2 ). 

The other alternative is simply that the Bondi accre- 
tion rate estimates at large radii are in error, with the 
true accretion rate being <5 10 _5 Mq yr -1 . We see several 
possible reasons for such an error: 

1. The estimates of the stellar mass loss parameters, M w 
and V w , may be in error. Note that a 25% uncertainty in 
V w translates to a factor of 2.5 uncertainty in the Bondi 
accretion rate. Together with a factor of few uncertainty 
in the stellar mass loss rates, M w , the Bondi accretion rate 
is probably uncertain by at least a factor of ~ 5. 

2. If the stars in the Galactic Center are not randomly 
distributed around Sgr A*, but instead have a large z- 
coordinate offset, the predicted accretion rate will be re- 
duced. We find that a systematic offset of <; 0.5 pc along 
the line of sight is needed to bring M a 1O _5 M yr" 1 . 

3. Perhaps Sgr A* has a strong outward wind which 
maintains a low-density "bubble" around the source (the 
winds from the stars in the small stellar cluster within 1" 
of Sgr A* - possibly 09 stars; Genzel et al. 1997 - might 
also contribute sufficient momentum to partially impede 
accretion of the Hel stars winds). If this bubble extends 
beyond the accretion radius Ra, then the accretion rate 
could be significantly reduced. The bubble is likely to be 
held up by mechanical pressure from a wind rather than 
radiative heating since the luminosity of Sgr A* is much 
too low for the latter to be efficient. 

We thank Rob Coker for useful conversations and the 
referee for useful comments. We acknowledge support 
from an NSF Graduate Research Fellowship (E.Q.) and 
NAG 5-2837 (R.N.). 
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